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Abstract 

Fitting Z^-pole data from ALEPH and SLD, and TPC data at a lower c.m.s. energy, we 
fix the boundary condition for NLO parton— s>hadron (hadron= vr^, K^, h^) fragmen- 
tation functions (FFs) at the low resolution scale of the radiative parton model of Gliick, 
Reya and Vogt (GRV). Perturbative LO^NLO stability is investigated. The empha- 
sis of the fit is on information on the fragmentation process for individual light (m, d, s) 
and heavy (c, b) quark flavours where we comment on the factorization scheme for heavy 
quarks in e~^e~ annihilations as compared to deep inelastic production. Inasmuch as the 
light quark input parameters are not yet completely pinned down by measurements we 
assume power laws to implement a physical hierarchy among the FFs respecting valence 
enhancement and strangeness suppression both of which are manifest from recent leading 
particle measurements. Through the second Mellin moments of the input functions we 
discuss the energy-momentum sum rule for massless FFs. We discuss our results in com- 
parison to previous fits and recent 3-jet measurements and formulate present uncertainties 
in our knowledge of the individual FFs. 



1 Introduction 



In this article we consider the production of identified hght hadrons in e~^e~ colhsions or 
more precisely the production of charged hadrons built from light up, down and strange 
valence quarks, i.e. dominantly pions, kaons and nucleons. Within perturbative QCD the 
light hadron production dynamics of e~^e~ annihilations inevitably comprise a nonpertur- 
bative long distance component from the nonperturbative hadronization of perturbatively 
produced partons. We will fix this latter nonperturbative component in the canonical 
QCD framework of parton^hadron fragmentation functions (FFs) at the low resolution 
input scale of the radiative GRV parton model In the naive parton model the FF Dp[z) 
has the interpretation of a probability density that some final state parton p hadronizes 
(or fragments) into a mean number D^{z)dz of hadrons h per dz where z is the fractional 
momentum which h receives from the parton. Field and Feynman have constructed |0| an 
early set of mesonic fragmentation functions footing solely on this intuitive probabilistic 
interpretation. Today's state of the art embeds the Dp{z) functions in the framework of 
QCD factorization theory ||^ including explicit operator definitions 0|. Fragmentation 
functions are the final state analogues of the initial state parton distribution functions 
(PDFs) and precisely as the PDFs do the FFs parametrize our ignorance of QCD bound 
state dynamics. Fragmentation functions into identified hadrons are therefore in their 
own right an interesting source of information on the hadronization process and they are 
a necessary ingredient to interpret present and future measurements of any one-particle- 
inclusive hard cross section within fundamental (perturbative) QCD theory as compared 
to Monte Carlo model approaches. Whereas PDFs are mainly determined from fully 
inclusive DIS the cleanest extraction of FFs is from e^e~ collisions. FFs pinned down 
in e^e~ can, by universality, then be applied to, e.g., the hadro- or leptoproduction of 
identified hadrons and compatibility with the transverse momentum spectrum produced 
in photon-proton collisions was found in |Q . From a combined use of both - fragmenta- 
tion and parton distribution functions - in one hadron inclusive deep inelastic scattering 
information on the initial state parton fiavour can be obtained from leading particle ef- 
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fects where a high energetic hadron inside a jet remembers the valence parton it has been 
produced off Parton information from semi-inclusive measurements based on an 

understanding of the fragmentation process also extends to polarized DIS where the spin 
flavour structure of the nucleon sea is still rather unknown [§,0. In this article we will 
restrict ourselves to the deflning e~^e~ process and leave other applications to future work 



The recent years have seen much effort p|,pJ]- [l^j to establish a similar technical skill 
for FFs as for existing sets of PDFs [1,1^, 2^. We will take a further step towards this 
goal and for this purpose we concentrate mainly on data which constrain the ffavour 
decomposition of fragmentation spectra; such data sets exist at the pole for charged 
TT^ and mesons from SLD ||2l|] and for the inclusive sum over charged hadrons from 



ALEPH ||Tj]. Lower energy counterparts from TPC will be included to correctly 
take QCD scaling violations - which have been re-established recently in |^,|13|,|l7| - into 
account. The flavour separation confronts us with the question of how to treat heavy 
flavours within e~^e~ production dynamics which we will answer in some detail. Especially 
the SLD data sets on flavour tagged fragmentation into vr^ and mesons contain new 
flavour information not included in previous flts to which we compare our results, 
thereby formulating the present uncertainty of our knowledge on the individual FFs. 



2 Parton Fragmentation in e+e Collisions Beyond 
the Leading Order 

The NLO framework^ for one-hadron-inclusive e^e~ annihilations is well known since 
long [p|, pn], P^, |2B|-p7| and we restrict ourselves to a brief theoretical introduction here 



in which we closely follow Ref. |11] in notation. 



^The formulae below include the LO framework in an obvious way by dropping subleading terms. 
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To be specific we will consider the reaction 

dz dz dz dz ' 

where 

z = 2Eh/Q = 2Ph ■ q/Q^ (2) 

is the energy^ Eh of the observed hadron scaled to the beam energy Q/2 = ^/i/2;0 with 
the positron/electron beam momentum Pe± = {Q/'^i 0, 0, ±(5/2) and q = + P^-. The 
r.h.s. of Eq. (|l[) distinguishes the contributions from transverse (T) and longitudinal (L) 
virtual bosons, where the polarization axis is in the direction of the momentum of the 
observed hadron h. Experimental data are commonly not presented for the absolute cross 
section in Eq. (|^) but for the normalized distribution 

1 AN^ I da^ ^ r. 

W- -A ' — ; as A^; ^ , 3 

Ntot Az atot dz 

where AA^'^ counts the registered /i-events per bin Az and Ntot denotes the inclusive sum 
of hadronic events. Within perturbative QCD and up to 0{al) the total hadronic cross 
section atot is given by 



O ( -g^J (4) 



where the parton model, i.e. 0(0;^), electroweak cross sections ctq for producing a qq 
pair are given in the Appendix. The power suppressed terms arise from perturba- 
tive quark mass effects [0{m?^/ Q"^)] or nonperturbative higher operator matrix elements 
[0[Kq(jj^/Q'^)]. We will discuss the quark mass effects in Section ^ where we give ample 
reasons for our choice of the factorization scheme. Higher operators are not considered in 
this article. 



^ In cases where experimental data are presented in the momentum scahng variable Zp = 2ph/Q we 
will transform the data to z using the relativistic energy momentum relation Ef^ = Ph~^ ^"h ^"^^ assuming 
pion production dominance (m/i — ttItt) whenever the hadron h is not specified. 

^ Though s = we will use both variables - s and - in the following in their role as the c.m.s. 
energy and perturbative hard scale, respectively. 
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QCD factorization theory predicts that 



dz 



E 



q — li, d, s, 
q — d, s, 
9 



dC 



2 



(5) 



+ 



A 



where the power n of the nonperturbative corrections ||28| cannot be determined from 
an operator product expansion based analysis, contrary to = 4 in Eq. (^) and to deep 
inelastic scattering where higher twists are known to be suppressed by a power n = 2. As 
anywhere the factorization and renormalization scale, respectively, and we will 

set them both equal to /i^ = /i^^^ = in the applications. The treatment of charm and 
bottom in the sum over quark flavours q = u,d, s, ... in will be discussed in Section |^. 
The coefficient functions Cp are given up to 0{al) in the MS scheme |TT],pB|,P^ by 



O P Crp 



27T 



27r 



CfcI 



Cp {CQ'^, iA^,r) 

Cp {CiQ'^i f^'F,R) 

CliCQ^f^ln) = ^4!^C,c^,(c,^)y:a^,{s) (6) 

and for antiquarks Cp p = Cp ^. The Cp^p and Ug can be found in the Appendix. The 
0{a1) contributions to the coefficient functions in are known |]29| but are of 
next-to-next-to-leading order and will therefore not be considered in this NLO analysis, 
with one exception to be discussed in Section 0. 



In the MS scheme the parton model expectation 



[\zzJ2 Dl{z) = l 

Jo u 



(7) 



that the entire parton energy is shared by the parton's fragmentation products is preserved 
under renormalization [D'l{z) /i^)] of the fragmentation functions by energy 



conservation 

1 . ^ da^ 



h 



where da^ is related to the via Eqs. (|T|), (|^). The renormahzed fragmentation functions 
obey massless AltareUi-Parisi-type renormahzation group equations 



91nQ2 

where the Pij have a perturbative expansion 



(9) 



P,(.,Q^) = ^ (.) + (^)' P«(.) + O . (10) 

The NLO pieces P-p of the timehke Pij in ([T0| ) differ from their spacehke counterparts. 
The Pij are imphcitly understood to represent the timelike sphtting functions in ||25|| P|. 

The evolution equations will be solved analytically in Mellin n-space as described, e.g. 
in [^^. We have included the n-space expressions |]2^,|2^ for the c^'?^ in Eq. in 



the Appendix. The timelike splitting functions ([T0|) up to two loop order have been 
transformed to n-space in [^]. To keep this article compact we do not reproduce these 
lengthy formulae nor do we repeat the solution of Eq. in Mellin space. Suffice to say 
here that some functional input forms for the D'^{z,fi'^) are required where we make the 
Ansatz 

Dl{z,f^l) = N,z-'^ {l-zr"^ (11) 

which we assume to hold for light partons (i = g,u,d,s and corresponding antiquarks) at 
the low input scale /Xq = 0.4 GeV^ of the recent NLO revision |jl| of the radiative parton 
model p0|^3|. (Accordingly, fi^ = 0.26 GeV^ |l[] will serve as the input scale for the 



LO fits to be discussed below.) As already mentioned, the treatment of heavy fiavours 
will be specified in the next Section. Along with the low input scale fi^ we also adopt 
the evaluation of af^^*^ used in |jl|, i.e. we numerically solve the renormahzation group 
equation 

^ aliQ') - aliQ') (12) 



dlnQ2 ' i67r 



*A misprint in the second Ref. of psj was corrected in |26|,P7| 



exactly by finding the root of 
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47r 



Pi. 



(13) 



PoasiQ^) 

witli /9o = 11 — 2//3 and /?i = 102 — 38//3 and wliere tlie number / of active flavours in 
the quark loop contributions to the beta function is 



/ 



(14) 



with |1|] rricfi^t = 1-4, 4.5, 175 GeV. Eq. (|T4|) guarantees the continuity of at the 
transition scales mc^i„t in the MS scheme up to two loops , if we furthermore adopt |1|] 
A3,4,5,6 = 299.4, 246, 167.7, 67.8 MeV in NLO (and A3,4,5,6 = 204, 175, 132, 66.5 MeV in 
LO). The continuity of at the transition scales Qq, where + 1, is guaranteed for 

Qq = m^^^ from the ratio of the renormalization constants Z3 in a renormalization 



scheme with / and / + 1 active flavours, respectively. Hence the continuity is not affected 
by the choice of solving the NLO renormalization group equation exactly as implied by 
Eq. (|T^) or by an analytical approximation up to some inverse power of InQ^, as, e.g., in 
Eq. (9.5a) of |^. The exact numerical solution of Eq. (^) is more appropriate |jl|] in the 
low < ml regime and will be used over the entire range. 



3 The Longitudinal Cross Section da^/dz 



We shall comment here briefly on the counting of perturbative orders for the longitudi- 
nal structure function on the r.h.s. of (|1|). Experimentally da\/dz is extracted by 
reweighting inclusive hadronic events according to their polar angle {6) distribution as 

(Pa'' 



dz 



dcosO 



dzd cos 6 



WLicose,v) 



where the projector Wl was introduced in [|ll| and v = \ cos 9 



(15) 

is set by the detector 



geometry, e.g. v = 0.94 for ALEPH 
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From the theoretical side we can read off Eqs. (D, (|) that dal/dz receives its leading; 
nonzero (finite and scheme-independent) contribution at 0{a\). When considering data 
on da\ obtained from reweighted events as in (|15D we will therefore include for our next- 
to-leading order analysis the 0{a1) contributions of Ref. We will, however, treat 
the 0{a\) coefficients C|^'^ in (|]) as subleading (NLO) contributions to the total da'^ = 
dal^ + da^ in Eq. (|I|) which receives its leading contribution from the da^ component to 

To 0{a^^ the convolutions in (|^) are most conveniently decomposed according to the 
flavour group as described in |2^. The corresponding coefficient functions for da^/dz are 



given in Eq. (17), (18) and (20) in the second Ref. of p9[. For the transformation to 



Mellin n space the nontrivial ones of the required identities can be found in [36- 



4 Treatment of Heavy Flavours 

Heavy quark initiated jets provide a substantial contribution to one-hadron-inclusive e+e~ 
annihilation spectra. In principle, the hard scattering production of heavy quarks has a 
well defined and unique perturbative expansion within QCD but the residual freedom 
of arranging the perturbation series at finite order leads into scheme ambiguities . 



The key question is whether the convergence of perturbative calculations improves if 
fixed order calculations are supplied with an additional resummation of quasi-coUinear 
logarithms [{as/^Tr) ln((5^/m^)]" to all orders n, where m is the heavy quark mass. The 
explicit all order resummation is equivalent to solving massless evolution equations from 
an input scale of 0{m). The boundary condition at the input scale is again calculable at 



fixed order and has been derived for e^e annihilations in |50] to 0{as)- 



Fixed NLO calculations seem at present to be reliable [jl],^ for the deep inelastic 



production of charm and bottom - dominantly in scattering events off wee gluons at low 
partonic c.m.s. energy where mass effects are most important [^. Contrarily, in e^e~ an- 
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nihilations each quark of a primary heavy quark pair receives the energy (5/2 3> m at the 
intermediate boson decay vertex. Such a highly relativistic quark will obviously not 'feel' 
its mass anymore and at the energies which we will consider, primary up- and down-type 
quarksQ are each produced in equal number, irrespective of the quark mass. We therefore 
expect heavy quarks in e~^e~ to behave essentially like massless partons with mass (quasi- 
) singularities to be resummed. In |^2|- ^] all order NLO collinear resummations were 
demonstrated to be required and adequate in order to describe the energy distribution of 
charm and bottom quarks over the partonic final state of high energy e~^e~ annihilations. 
In particular, all order massless resummations for charm quarks were proven to be nec- 
essary to describe the amount of secondary charm quarks from (7 — > cc splittings in the 
parton showering, visible in a rise of the cross section at lower z as observed by OPAL 



and ALEPH |^ and as opposed to the fixed Ola"^) order contribution qqg — ^ qqcc 



47| which turns out to be far too small p3|,|4^. Bottom production has an evolution 
length which is shorter by the amount of AlnQ^ = ln(m^/m^) which suffices to suppress 
secondary bottom pairs as experimentally observed in a flat z ^ B-spectrum. 

Hence, we will include q^ ^ h FFs mixing under evolution with their light parton 
analogues, i.e. describing all long distance (collinear parton showering and hadronization) 
effects of the fragmentation process. Heavy flavours are treated above their respective 
MS 'thresholds' Qq = rric^b as active flavours in the evolution of the D^g^z,^!^) in Eq. (H). 
We adopt the functional form of Eq. (|TlD also for D^^, i.e. 



Dl,^,{z,^^^) = N,z^U^-^f , (16) 



to hold at /if = mf along Eq. ([T^) which guarantees the continuity of as at fif in NLO. 
The D^^^fj are then c?iscontinuously 'switched on' in the evolution at rrif but enter the 
cross section in Eq. (|l]) only above the partonic threshold > 4m? .| It is certainly 
questionable ||13[ to fit pure QCD fragmentation functions for heavy quarks to inclusive 



^ We do, of course, not consider the superheavy top quark here. 

^ The partonic threshold 4mf - or similar low resonance threshold scales [|j - may therefore be 
considered a natural choice for fif. We follow existing PDF sets which avoid the discontinuities in as 
induced by fj,i ^ rrii. 
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or heavy quark tagged data because charm and bottom jets are highly contaminated by 
weak decay channels. Adopting such a procedure anyway is a necessity as long as data 
corrected for weak decays do not exist. 



5 Fitting Procedure 



Early model attempts to generate fragmentation functions entirely by QCD dynamics 
from a delta peak input Dp{z) oc 5(1 — z) are outruled by modern high statistics data 
|l3,|2l],^ which require smooth input functions (|lTl), (p^Bf ) to be fitted to experiment. 
Recently, fragmentation functions have been fitted to e~^e~ production data using either 
free fit Ansatze or Ansatze strongly constrained by SU(3)f symmetry |T^. We will take 



an intermediate path and constrain free fit Ansatze by making power law assumptions 
about valence enhancement and strangeness suppression within the radiative parton model 
of Refs. [P], P(]|-|55[| . The model is originally tailored for parton distribution functions 
to perturbatively generate the high population of quarks and gluons in hadrons at low 
Bjorken x. It should be clear from the scratch that its predictivity cannot be transfered 
to FFs where a well defined low z analogue of the deep inelastic lovj x regime is missing 

Our aim is to pin down fragmentation functions D'^z, ij?) for (anti-)quarks and gluons 
hadronizing into charged pions {h = 7r+,7r^ = 7r+'^), charged kaons {h = K^^K^ = 
K^~) and the inclusive sum over charged hadrons [ Ylhi^^ + ^ ) = Xl^^lsQ basically 
from high statistics data at the peak measured by ALEPH [T^ at LEP h^) and by 
SLD {tv^, K^) at SLAG. Besides their high statistical accuracy, these data sets have 
the advantage of furnishing along with their fully inclusive measurements also flavour 
enriched events which can be used to fix the flavour structure of the D^{z,^'^) to some 
extent. Both data sets distinguish between light quark (m, (i, s), charm, and bottom events 
where the quark flavour refers to the primary Z°-boson decay vertex. Since pure flavour 
^ Similarly, tt^ will denote 7r+ + tt^ to be distinguished from 7r+'^ = 7r+,7r^; for kaons accordingly. 
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separated sets cannot be obtained directly |T7|,pl|, Monte Carlo simulations are required to 
estimate the flavour composition of the tagged data sets. From these Monte Carlo studies 
the flavour enriched event samples have been unfolded at a given systematic uncertainty 
to pure flavour sets for the SLD data whereas ALEPH quotes percentages for each flavour 
contribution to any of its flavour enriched samples. In the latter case some uncertainty will 
reside in the translation of the Monte Carlo studies to the perturbative calculations we 
are performing here which can unfortunately not be quantified since no systematic errors 
are quoted []T^ for the percentages. Anyway, we will reweight the electroweak couplings 
in Table in the Appendix to reproduce the ALEPH flavour composition as quoted in 



17|. In order to have our EEs respect QCD scaling violations properly we include lower 



scale {^/s = 29 GeV) TPC data p2[ in our fits which also furnish flavour information 
(unfolded to pure {u, d,s}, c and b event samples) for ^ and for tt^ but only inclusive 
measurements for K^. Needless to say, the fits will be dominated statistically by the 
Z^-po\e measurements of ALEPH and SLD. 

Our Ansatze for the fragmentation functions will be the ones of Eqs. ([TT|) , ([T6|) for 
light parton and heavy quark fragmentation, respectively, where we assume the following 
symmetries and hierarchy 

Df^h-^jjh-M. h = ^^x (17) 
Df = < Df = Df (18) 
Dr = dS < nr < dT , (19) 

where Eq. ( p!7[ ) respects charge conjugation and Eqs. (p!8|), ([T9|) should hold from the 
valence structure of pions and kaons and strangeness suppression. The equality in Eq. 
(p!8|) seems to be confirmed by 'leading particle' measurements |^0| which also indicate the 
second mequality in Eq. (|l^) at large z from the suppression of secondary ss formation 
which is required to form a from a u but not from an s quark. Beyond strangeness 
suppression, in 7r+'^ or production we expect differences from favoured valence- 

type (e.g. u 7r+) and unfavoured sea-type (e.g. s 7r+) fragmentation channels. 
Nevertheless, we assume a universal z ^ behaviour, determined by the input parameters 
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Ni and in (|TT]), This assumption is guided by the idea that the valence enhancement 
of, say, M ^ 7r+ fragmentation should manifest itself mainly as a 'leading particle effect' 
parametrized by (3i in (|ll|) as 2; ^ 1. We will assume |^ 



which suppresses ss formation as well as sea-type fragmentation as z ^ 1. A linear 
suppression factor DJ^/D^^ = {1 — z) for sea-type fragmentation is compatible with 
semi-inclusive deep inelastic measurements ||^,|^ in the range 0.1 < z < 0.8. These 
measurements seem to prefer a factor ~ (c — 2;) with c ~ 0.25 for z > 0.8. The very large 
z behaviour of the D^{z) is, however, not very well constrained by e^e~ measurements 



and from the theoretical side soft gluon resummations may become necessary |T8|. Also, 
the universality, i.e. process-independence of the FFs in e^e~ annihilations and semi- 
inclusive DIS|[T^, respectively, is not yet well settled at large z We therefore keep 
D'^^ / D'^ = {1 — z) for the time being. Our assumptions in Eq. pO| ) are compatible 
with DJ^ {n = 2) ~ 0.6 D'^^ {n = 2) where the second moments are defined below in Eq. 
(p^) and where their ratio can be estimated from EMC one-pion- inclusive data 0. These 
data have been analyzed in a parton model context and the extracted scale-independent 
FFs of can therefore not be compared to our fit in more detail. The corresponding 
one-kaon-inclusive measurements in |^ seem to prefer an even stronger suppression of 
sea-type fragmentation in kaon formation. For the time being, we do, however, assume 
the suppression of sea-type channels to be a universal phenomenon modeled by one extra 
power in {1 — z) for the input FFs. This assumption, as any of our above assumptions on 
the light flavour structure of annihilation data, must obviously be expected to be violated 
to some extent if dedicated measurements will become possible. 

If we furthermore assume for simplicity 

= D^- (21) 

to hold in the heavy quark sector^, we are left with the following set of independent input 

^ We would expect Eq. (pi|) to hold exactly if heavy quark jets were not contaminated by weak decays. 
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parametrizations for pions and kaons: 





= Nfz<^ 


(1 


- ^)^"^ 




= Nfz<^ 


[1 






= Nfz^'f 


[1 


-zff 




= Nfz'^-^ 


[1 


- zT^^ 




= Nfz<^ 


(1 


- z)f'f 



Dr 






= iVf+z-" 


(1 


— z 


fV 




{z 




= N^\"^^ 


(1 


— z 






{z 




= N^z-^" 


(1 


— z 










= Nf^'^?^ 


(1 


— z 




c,c \ 




ml) 


= iVf^""^ 


(1 


— z 




^b,b \ 


z, 


ml) 


= 


(1 


— z 





(23) 



In practice, we will express the normalizations through the physically more interesting 
contributions to the energy integral in Eq. (0), given by the second Mellin moment 



Df(n = 2,/x2)^ [\zzDl{z,f,' 
Jo 



) 



(24) 



which would be an exact invariant under evolution in fi"^ only if D^{n = 2,yU^) were the 
same for all i - which is not the case. Anyway, guided by the idea of having an intermediate 
gluon in the parton cascade connecting a valence- and a sea-type quark we will reduce 



the parameter space {N^ 



TT+-,K + ' 



, j3j ' ' } somewhat by demanding 



Df{n = 2,f,l) 
Df{n = 2,f,l) 



Dfin = 2,f^l)+Dfin = 2,fil) 
Dr{n = 2,f,l) + Dr{n = 2,^^l) 



(25) 



which reduces the strength of the evolution of the D'l{n = 2,yU^) considerably and does 

not influence the quality of the fits. 

® We do not include in the second average of Eq. (|2^) because gluon fragmentation into a K'^ 

must be strangeness-suppressed. 
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Samples of inclusive charged hadrons ^ h"^ are dominated by charged pions and adding 
a kaon background is a reasonable approximation for most applications: 



A - A • (26) 

Eq. (^) is, however, not adequate to compare with high statistics measurement from the 
Z°-pole |r^, which are sensitive at the (9(l%)-level. We will fit a small residue from 



higher mesons and baryons along 



= , (27) 



using again simple Ansatze as in (|TT]), (|16]) where we distinguish between light flavours 
{D'^''- = D^J'''- = Df""-), heavy charm and bottom quarks {Dl''''- , Dl''''-) and gluons (-D^'''' ) 
and where Dl"''- > serves as a consistency check. The residual functions amount to 
a sufficiently small correction not to constrain them any further. For the same reason 
the approximation that the hadronic residue is half positively, half negatively charged for 
any parton, i.e. Df^^' = D^^' + Df''- /2 can be safely used for the inclusive sum 

of only positively (negatively) charged hadrons. Altogether, the Ansatze (p2D-(p3D, (p7|) 
with the constraints ( p5|) will be fed into the evolution equation (|^) to evaluate the cross 
section (|l|) using the factorization theorem (|^). The x^-minimization algorithm MINUIT 
5^ will be used to obtain best possible agreement of the outcome with experimental 



data [p!7|, pT|, p2| over the range 0.05 < z < 0.8. Towards low z the TPC data for charged 
pion production ||22| at ^/s = 29 GeV lie above the corresponding SLD data [^] at 

^ = Mz (c?f^^^^L<oi > ^^^^°L<oi) ^^i^^ contradicts the established QCD 
evolution {da^'"^ /ds\_^^^ ^ > 0) predicting that da increases with ^/s at low z. Since the 



SLD measurements are compatible with the LEP data of e.g. ALEPH [jT^ we exclude 
TPC charged pion data below 2; < 0.1. Data on the longitudinal structure function da^^^ 
available from ALEPH will be considered for consistency but are not weighted in 



the fit because we found no signiflcant statistical impact of these data on the flt results. 
Although one might hope daL to constrain the gluon fragmentation function because the 
gluon enters at 'leading' 0{a].) a closer NLO inspection reveals that dai is dominated 
by quark fragmentation over most of the z range. 
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0.264 




A'-o) 


Mfz-^^-^'^l - 


^) 1.949 
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jyres. ^—0.605 ^j^ 
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Table 1: Numerical values for the NLO fit parameters in Eqs. (p^)-p3|), (|27D. The input 
scales are /iQ,m^,m^ = 0.4, 1.96, 20.25GeV^ |jl[. The Normalizations N^^ are determined 
by the second moments of the input functions. These moments evolve rather mildly 
(typically less than 10% up to M|) to higher scales due to the constraint (p5|). 



6 Results and Discussion 



The quality of the fit can be inferred from Figs. |1], |^, and 0. All non-6 event samples 
can be well described within errors, even the (non-6-tagged) ALEPH data where the 
error is dominated by a normalization uncertainty of only 1%. As already noted in 
|]T^,|S3| the b quark fragmentation spectrum can not be described perfectly well using the 
simple functional form (|l^) as input to the evolution. The b sets are only reproducible 
at some 5% accuracy. With the inherent uncertainty due to weak decays we consider this 
a reasonable precision and do not investigate more involved functional Ansatze for the 
-D^. Note that the longitudinal cross section (lowest curve in Fig. |I|), calculated along 
29| , has not been included in the fit. We list in Table |I] the input functions introduced 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




Figure 1: ALEPH [|T3 ^ /i^ inclusive particle spectra, measured at the Z° pole, and the 
corresponding fit results. Details to the individual data samples and curves are given in 
the text. The 'longitudinal' set has not been included in the fit. 
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Figure 2: SLD tt^ and inclusive particle spectra, measured at the Z° pole, and 
the corresponding fit results. Details to the individual data samples are given in the text. 



in Eqs. (p^ ) -(p3| ) , (|27|) along with their contributions to the energy integral As a 
representative we show the set Df'^ of input fragmentation functions into (positively) 
charged pions in Fig. ^. Evolution effects to higher scales can be inferred for the low 
input scale (/ig = 0.4 GeV^ functions -DJ^^ ^ fro^i Fig. ^. They are quite dramatic 
for the gluon FF but less pronounced for the light flavours. For the higher scale (^ct) 
input functions -0^6) shown in this Figure, the evolution effects are of course still 
weaker as for -DJ^^. The peculiar shape of the input Dg^{z,fil) should be traced back 
to starting the evolution at the low input scale of ^ thereby maximizing the driving 



16 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



10- 



10 

10" 
10 
1 




i_ inclusive: do/dz 
I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




uds events :(|3o)' do/dz 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i TT7 f I I I I 




0.35 2 0-4 



0.45 



Figure 3: Yl ^^"^ inclusive particle spectra and flavour separated events over the 
range < z < 0.9 (inclusive) and < z < 0.45 {{u, d, s}, c, and b events; /5 = Pn/E^^ — 1 
for not too small z) as measured at ^/s = 29 GeV by TPC |22]. The corresponding curves 



are the fit results. Details to the individual data samples and curves are given in the 
text. Also included in the fit were inclusive data by TPC which are, however, not 
accompanied by flavour tagged samples and are therefore not shown in the Figure for 
clearness. 



force of Dg . From pion dominance we would assume a value of Df^' {n = 2) < 2/3 
and it is interesting that this expectation is roughly confirmed. Also, estimating the 
contributions from neutral hadrons from SU(2)f symmetry results in values for the r.h.s. 
of Eq. within [0.9; 1.0] which confirms energy conservation (H) except for Df^ which 
violates to about ~ 10%. However, as mentioned above, a pure QCD fit to bottom 
fragmentation is questionable due to weak decay channels which mostly contaminate b 
jets. Furthermore, a physical interpretation of the second moment integrals (^) over the 
range ze[0] 1] is delicate and may be misleading in general due to a sizable contribution 
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Figure 4: The input fragmentation functions D^"^ as given in Table ^ at their respective 
input scales. 



from the perturbatively unstable very low z region, as we will discuss in more detail in 
Section lOI. 



We have also evolved our ^ fit to a c.m.s. energy of yfs = 161 GeV and compared 
the evolution to measurements by OPAL [5^| in Fig. |^. The agreement is convincing albeit 
not too surprising because scaling violations from Mz — > 161 GeV are rather moderate. 
In Figs. 1^ and |^ we compare our fitted fragmentation functions to a previous NLO fit by 
Binnewies, Kniehl and Kramer (BKK) 0. These authors confirm QCD scaling violations 
within a wider range of data, albeit not with the flavour information of Ref. on the 
individual vr^ and spectra. Instead, the flavour separation in |^ was done by fitting 
to the flavour tagged ALEPH [T^ ^ data and assuming da^'' = da^^+^^ + / where 
/ is a small residual from (anti-)protons as measured in |5^. Besides discrepancies in the 
barely constrained D^, our procedure of decoupled fits to Yl^^^ ^^"^ data yields 
quite different results for the individual flavour fragmentation functions into vr^ and K"^. 
Note that the differences of our fit to BKK decrease stepwise the more fiavour-inclusive 
sums of FFs are considered. For the 'democratic' FF D^^^^^^^^^^^ the difference shrinks 
to at most 4% within 0.1 < z < 0.8. In any case, the light flavour {uds) structure is, at 
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Figure 5: The input for light {uds) quarks and gluons of Fig. ^ evolved upward to = 
Ml 

present, arbitrary to some extent in our fit as well as the one of Ref. P| and differences 
between the two fits of the individual DfJ}^^ ^ estimate the present uncertainty on these 
functions. 



6.1 3- Jet Measurements and Gluon Fragmentation; Leading Par- 
ticle Effects 



Complementary to our considerations of one-particle inclusive e^e" annihilation data, 
individual quark and gluon fragmentation functions can also be derived from 3-jet topolo- 
gies which can at the leading perturbative order be attributed to partonic qqg 



configurations and where an experimental fragmentation function can be defined as 

1 A jvfh rp 

i.e. as the spectrum in energy Eh scaled to the jet energy Ef^ of some hadron species h 
distributed inside a jet initiated by a parton i. In general, QCD scale dependence from 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 z 0.8 

Figure 6: OPAL inclusive particle spectrum measured at a c.m.s. energy of ^/s = 

161 GeV. The curve is the fit of Fig. ^ evolved upward to that energy. 

leading logarithms can be understood non-covariantly from the transverse phase space 
volume available for collinear parton emission {splitting process) 



dpi 

pI 



(29) 



and a phase space boundary {p^^^^ = 0{s) leads to the common choice /x = = Q 
for an (one-particle-) inclusive phase space. A different p™*^^ is, however, induced in a 
three jet topology from the requirement that parton emission proceeds into a cone the 
geometry of which is defined by the hadrons which are grouped together as the jet. The 



analysis in [58| has demonstrated the fragmentation function defined in Eq. (28) to un- 
dergo scaling violations compatible with LO QCD evolution in the jet topology scale 
/i = E^^^ sin{6/2) where 6 is the angle to the nearest jet in a 3-jet event. For the time 
being, these considerations have to be restricted to a LO treatment, since NLO corrections 
(stemming e.g. from a kinematical configuration where the quark of the qqg triple emits 
a high pj_ gluon into the jet clustered around the anti-quark) have to our knowledge not 



been formulated yet and we have therefore refrained from including the data of [56-5S 



in our fits. The gluon fragmentation function is, however, barely constrained from 
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Figure 7: Ratios of the individual fragmentation functions obtained from this NLO fit to 
their analogues in Ref. 

one-particle-inclusive final state measurements in e~^e~ annihilation because it enters the 
cross section in Eqs. @, (|^) only in subleading order 0{al) and where the leading part 
is factorized into the evolution where mixes with the quark singlet fragmenta- 
tion function D^. Therefore, keeping the vagueness of this comparison in mind we 
use the 3-jet data of Ref. to compare our fitted D^^^ [z, with, i.e. we compare 
D^^^ (^z = XE,fi'^ = {Ef^f ^ (40GeV)2j with the measured ® 
time being, with the caveat that this cannot give us more than an idea of the compati- 
bility of our NLO gluon FF with 3-jet measurements. For illustration we also include the 
gluon fragmentation function of [using the approximation (|26|) ] and an independent 



experimental LO determination from DELPHI |£8[ in our comparison. Hence, the dis- 
crepancy between the data |^ and the LO QCD fit to an independent measurement 
estimates the accuracy to which the 3-jet data can be identified with a QCD gluonic 
fragmentation function. From Fig. |^ we judge that the result seems promising and that 
a refinement of the theoretical framework would probably contribute to removing the 
existing ambiguities in the gluon fragmentation function. 



Similar theoretical limitations as outlined above for 3-jet measurements also prevent 



21 



3.5 
3 

2.5 
2 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



u, s 
d 
c 
b 



= 100 GeV^ 



d'^(z,q'): 
this fit / BKK 




Figure 8: Same as Fig. ^ but for the fragmentation functions into charged Kaons K"^. 



at present a detailed pQCD analysis of leading particle effects [£0[ in e^e~ annihilations 
which are based experimentally on phase space restrictions which do not fully match a 
one-particle inclusive QCD framework. Qualitatively, the results of the leading particle 
measurements in are accounted for by Eq. (pO]). 



6.2 Perturbative Stability and Energy Sum Rule 

The full NLO framework of Section is expected to yield more reliable and less scale- 
sensitive results compared to a LO truncation of the perturbation series where the known 
0{as) terms of the coefficient functions (^) are neglected as well as the (3i contribution 
to the running of and where the omission of the P^^^ parts of the splitting functions 
reduces the evolution (j^) to summing only the most dominantly leading logs {asjln In Q^)"" 
for all n. Still, an accompanying LO fit is - besides future effective LO applications - 
valuable as a test of the perturbative LO^NLO stability which is a delicate requirement 
for perturbative (QCD) approaches to strong interaction phenomena - especially if the 
perturbative QCD dynamics is supposed to set in at the rather low input scale in ([TT|) 
taken from For infrared unsaie quantities - such as the one-particle- inclusive spectra 
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Figure 9: The gluon fragmentation function as defined experimentally in (pHj) and mea- 
sured by OPAL compared to the NLO fits of our analysis as well as the one in 



An experimental LO determination from DELPHI |5q| is also included for comparison. 



considered here - the nonperturbative parameters have to be redefined at each perturbative 
order since they replace new infrared sensitive terms of the factorized perturbation series 
at any order. We have accordingly performed an accompanying LO fit following exactly 
the same lines as the NLO fit described above and resulting in Table 0. In order not to 
overload the logarithmic-scale figures of the preceding sections with narrow pairs (LO, 
NLO) of lines and instead of repeating the details for the LO fitting procedure we rather 
concentrate in this separate Section on two LO^NLO issues worth mentioning. As a 
representative we will consider our sets of pionic fragmentation functions. In Fig. |10| we 



show ii'-f actors for da in Eq. (11]),/^ = c^ctnlo/ '^^lo' ^^^^ each flavour contribution 



Df + Df 



T+L 



(30) 



where C^.r = Cf 



T+L = ^T+L X ko/ Eg' ^0 ] is the Cf,^^ of Eq. 



I 'per flavour' and where 
denotes a convolution integral as in Eq. (^. Note that the 0{a].) terms of Eq. 
entering Eq. (|30| ) in NLO are neglected in LO where only the ~ 5(1 



z) term of 



contributes. The vertical lines in Fig. |10| indicate the range in z where the fit is rather 
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^^.+ ^-0.923 
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^) 1.976 


0.244 




A'-o) 


^5.271 ( 


1 - 


^^8.235 


0.311 






^^.+ ^-0.818 


(1- 


^)3.461 


0.241 






^^.+ ^-1.072 


(1- 


^^6.695 


0.264 


or 




^0.617 


;i - 


^)0.744 


0.213 


u 






;i - 


^) 1.744 


0.085 








;i - 


^)2.744 


0.044 






NK+ ^8.132 


1 - 


^)5.776 


0.064 


D^: 






;i- 


^)6.171 


0.085 


^b,b 




NK+ ^0.191 


;i- 


^^8.934 


0.062 


J~\TCS. 
1 




jyres. ^0.938 


(1- 


- ^) 7-734 


0.146 


jyccs. 

a 




jyres. ^6.150 


(1- 


_ ^^5.379 


0.003 






jyres. ^—0.636 


(1- 


- ^)2-486 


0.084 






jyres. ^—0.736 


(1- 




0.137 



Table 2: Input parameters as in Table |T] but for our LO fit where jUq = 0.26GeV^ [|T 



well determined, where z > 0.05 excludes the perturbatively unstable |jT6| low z region 
and where for z > 0.7 statistics drop rather low (only one data point for z > 0.7). 
Hence, the spread of the curves towards large z should not be taken as a perturbative 
instablitiy but attributed mainly to the decreasing experimental statistics which do not 
determine the FFs very well 



in the z ^ 1 hard fragmentation limit and where poorly 
defined minima may easily fake perturbative instability in decoupled LO^NLO fits. 
It is therefore reassuring to observe that the fiavour-inclusive da has a ii"-factor closer to 
one than the individual fiavour-contributions where the latter have poorer experimental 
statistics and larger systematic errors than the former. On the other hand, the unstable 
z behaviour is unaffected by the choice of input parameters and can be traced 
back to including/omitting the NLO pieces Pj^p of the splitting functions in the NLO/LO 
evolution. Indeed, our choice z^i^ = 0.05 appears already to be on the edge of perturbative 
reliability. The far steeper LO evolution has profound impact on the second moments in 
Table |^ which violate the energy sum rule (|]). We demonstrate this point in Fig. |lT 
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0.7 z 0.8 



Figure 10: K-factors for individual flavours and inclusive charged pion production. The 
vertical lines indicate the range where the fit is rather well defined where z^i^ is set by 
perturbative stability and ^max by experimental statistics. 



where we show the truncated moment 

, dzzVl\z,Q' 



(31) 



in LO and NLO. It can be seen that the full second moment receives in LO a sizable 
contribution from the unstable z < 0.05 region where the evolution generates a too steep 
spectrum resulting in large second moments violating (^). On the other hand, in NLO 
the second moment, truncated at Zjam = 0.05, is rather close to its wntruncated value 
leading to NLO moments that compare rather well with (^. One must note, however, 
that the choice z^i^ = 0.05 is rather arbitrary and that the region 2;e[0.0; 0.05] does not 
contribute much in NLO partly because (z,Q'^) turns unphysically negative at very 
small z due to the timelike NLO evolution [|T6|. This observation renders the energy sum 
rule (§) a delicate concept for perturbative QCD FFs and we believe it should not be 
considered within this theoretical framework unless the z — > behaviour of FFs is under 
better control. 
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Figure 11: The truncated second moment of the u-quark contribution to charged pion 
production as defined in Eqs. (|30|), (|3T|). Shown are results for our NLO and LO fit. 

7 Conclusions 



Within the radiative parton model [jT],^O|-0 we have fitted parton fragmentation func- 
tions to identified hadron (vr^, K^) and inclusive charged particle spectra measured at the 
pole |T7|,|21|. Scaling violations were properly taken into account by simultaneously 
fitting lower energy TPC data. Special attention was paid to the fiavour structure of the 
parton FFs where a collinearly resummed RGE formalism was argued to be adequate for 
the treatment of heavy quark contributions. In the light quark sector we made power 
law assumptions which qualitatively establish a physical hierarchy among the FFs which 
is guided by the ideas of valence enhancement and strangeness suppression and which is 
observed in leading particle measurements. It would be desirable to interpret these mea- 
surements - as well as gluon-jet measurements from 3-jet topologies - in more quantitative 
detail within the language of NLO QCD fragmentation functions. A corresponding frame- 
work has, however, to our knowledge not been developed yet. Despite the high precision 
data available, the individual Df^^'^ ^ '^^^ are therefore still rather uncertain. To estimate 



26 



the present theoretical uncertainty we compared our fit to the one of ^ and found siz- 
able deviations. An inclusive sum over the distinct flavours is, however, rather reliably 
determined for not too large z. The missing experimental information at large z may fake 
perturbative instability in independent LO<->NLO flts manifesting, however, just the fact 
that the FFs are still unknown in the z ^ 1 hard fragmentation limit. We also considered 
the contributions of the individual FFs to the energy integral in (^) and found reasonable 
values in NLO while the steeper LO evolution requires too large moments. The instability 
of the timelike evolution at low z makes, however, energy conservation a concept which 
may hardly be useful within perturbative QCD even for NLO FFs. 

A FORTRAN package of the LO and NLO functions will be obtainable upon request. 
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Appendix: 

MS Coefficient Functions for One Hadron Inclusive 
e^e~ Annihilation 



The 0{al) coefficient functions for One Hadron Inclusive e^e Annihilation as introduced 
in Section |I| read [p!l|, p^, [2^ : 



:i + c' 



[Ml -01 


3 


1 


[ i-C J 


+ ~2 





1 + 
1-C 



2l±il^rin(l-C) + 21nC + ln|)-4-^-^ 



c 



c 



4 



1-C 

c 



(32) 



The parton model, i.e. C(a°) electroweak cross section cTq for producing a qq pair in e~^e 
annihilation are given by 

4™^ 



[el + 2e,v,Vfp^{s) + {vl + a^)(t;J + aj)p2(s)] 



(33) 



with the QED fine structure constant a and where 

1 s(M| - s) 



P2K.S) 



2r2 ' 



4 sin^ cos2 9w (M| - s)2 + M|r| 



4 sin^ cos2 ) (M| - s)^ + M|r 



2r2 ' 



(34) 



and the electric charges Cj and electroweak vector (t>j) and axial (a^) couplings are listed 
in Table El according to their standard model values: 



Tg^j - 2ei sin^ 61,4/ 



- 3,i ) 



(35) 



where T is the weak isospin and 9w the Weinberg angle. 
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particle i 




Vi 


ai 


e~ 


-1 


-i + 2sin2 Ow 


1 

2 


up-type quark 


^3 


+i — 1 sin^ 6w 




down-type quark 


1 
3 


— i + 1 sin^ 6w 


1 

2 



Table 3: Electromagnetic and electroweak couplings entering Eq. 



The Mellin transforms of the c^^ in Eq. (|32D read 



n 



' 2 



5S2(ri) + S2(n) + SiH 



3 1 9 

2n + l n(n + l] 



1 



rafra + 11 



2Si(n) 



In ^ 



(n + l) 



n 



' 2 



-SiH 



+ n + 2 



+ 



(n — l)n(?2 + l) (n — 1)^ 



n 



' 2 
^ ' 2 



where the sums 



+ — TV^^~ 

n[n^ — 1) jjLp 

1 

n 



{n — l)n 



Skin) 



n ^ 



have to be analytically continued pll,p6| to the complex n plane 



SiH = 7f + V(^ + 1), 7£ = 0.577216 



S2H = C(2)-^'(n + l), C(2) 



(36) 



(37) 



(3^ 



with the help of logarithmic derivatives of the F function il)'^^\n) = d^'^^^^ \nT (n) / dn''^^ . 
Analogously to the 0{al) contribution Fl to the deep inelastic electron nucleon cross 
section, the 0{al) contribution from longitudinally polarized virtual bosons to the frag- 
mentation spectrum in e+e" annihilations is scheme independent and finite. The c^'^ do 
therefore not depend on the factorization scale contrary to the c^^ which are infrared 
safe only after factorization of the coUinear singularities. 
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